NIKHEF 04-004 
DESY 04-060 
SFB/CPP-04-12 
April 2004 



hep-ph/0404TTT 



The Three-Loop SpHtting Functions in QCD: 
o : The Singlet Case 

u _ 

; A. Vogt"^, S. Moch^ and J.A.M. Vermaseren'^ 



"NIKHEF Theory Group 
Kruislaan 409, 1098 S J Amsterdam, The Netherlands 



^Deutsches Elektronensynchrotron DESY 



O 

^ ' Platanenallee 6, D-15735 Zeuthen, Germany 

> 

X ■ Abstract 

We compute the next-to-next-to-leading order (NNLO) contributions to the splitting functions 
governing the evolution of the unpolarized flavour-singlet parton densities in perturbative QCD. 
The exact expressions are presented in both Mellin-A'' and Bjorken-x space. We also provide accu- 
rate parametrizations for practical applications. Our results agree with all partial results available 
in the literature. As in the non-singlet case, the correct leading logarithmic predictions for small 
momentum fractions x do not provide good estimates of the respective complete splitting func- 
tions. We investigate the size of the corrections and the stability of the NNLO evolution under 
variation of the renormalization scale. The perturbative expansion appears to converge rapidly at 
X > 10^-^. Relatively large third-order corrections are found at smaller values of x. 



1 Introduction 



Parton distributions form indispensable ingredients for the analysis of all hard-scattering processes 
involving initial-state hadrons. The scale-dependence (evolution) of these distributions can be de- 
rived from first principles in terms of an expansion in powers of the strong coupling constant ttg. 
The corresponding nth-order coefficients governing the evolution are referred to as the n-loop 
anomalous dimensions or splitting functions. Parton distributions evolved by including the terms 
up to order a"^^ in this expansion constitute, together with the corresponding results for the par- 
tonic cross sections for the observable under consideration, the N°LO (leading-order, next-to- 
leading-order, next-to-next-to-leading-order, etc.) approximation of perturbative QCD. 

Presently the next-to-leading order is the standard approximation for most important processes. 
The corresponding one- and two-loop splitting functions have been known for a long time El El 
m mm 1^1 [71 [HI muni nil- The NNLO corrections need to be included, however, in order to arrive 
at quantitatively reliable predictions for hard processes at present and future high-energy colliders. 
These corrections are so far known only for structure functions in deep-inelastic scattering (DIS) 
|[T2l[T3[[T4[[T5ll and for Drell-Yan lepton-pair and gauge-boson production in proton-(anti-)proton 
collisions |[T6[ [T7[ [T8[ [T9ll and the related cross sections for Higgs production in the heavy-top- 
quark approximation |[T7ll20l 121112211 . Work on NNLO cross sections for jet production is under 
way and expected to yield results in the near future, see Ref. ll23ll and references therein. 

For the three-loop splitting functions, on the other hand, only partial results had been computed 
until very recently, especially the lowest six/seven (even or odd) integer-N Mellin moments ll^E51 
1^ and the leading (Inx) jx small-x terms of three of the four singlet splitting functions liTTll^ . 
The results of Refs. JHEIIEEI have been employed - directly 1122111011111113 and indirectly lEH 
|33|| via Bjorken x-space approximations constructed in Refs. 0511^1771 from them and the small-jc 
constraints ll27ll28l - to improve the analysis of DIS data and hadron-coUider predictions. This 
information is however not sufficient for quantitative predictions at small values of x. 

We have recently published the non-singlet part of the unpolarized three-loop splitting func- 
tions ll38ll . In the present article we compute the corresponding singlet quantities. The article is 
organized as follows: In section 2 we set up our notations and very briefly discuss the method of 
our calculation. The Mellin-A^ space results are written down in section 3. The (In A'') jN subleading 
large-A'' term of the three-loop gluon-gluon splitting function is found to be related to the leading 
InA'^ contribution at second order, in complete analogy to the relation found for the non-singlet 
quark-quark case. In section 4 we present the exact results as well as compact parametrizations for 
the jc-space splitting functions and study their behaviour at small x. We demonstrate that neither 
do the (Injc) jx terms dominate the splitting functions at experimentally relevant values of x, nor 
do even all Xjx terms dominate the Mellin convolutions by which the splitting functions enter the 
evolution equations. The numerical implications of our results for the scale dependence of the 
singlet-quark and gluon distributions are illustrated in section 5. As in the non-singlet case the per- 
turbation series converges rapidly for x > 10^-^, while relatively large corrections occur for smaller 
momentum fractions. Finally we briefly summarize our findings in section 6. 
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2 Notations and method 



We start by setting up our notations for the singlet parton densities and the splitting functions 
governing their evolution. The singlet quark distribution of a hadron is given by 

'V 

1=1 

Here qj{x,ijj) and qj{x,ijj) represent the respective number distributions of quarks and antiquarks 
in the fractional hadron momentum ;c. The corresponding gluon distribution is denoted by g{x,ijj). 
The subscript i indicates the flavour of the (anti-) quarks, and n^: stands for the number of effectively 
massless flavours. Finally /Uf represents the factorization scale. For the time being we do not need 
to introduce a renormalization scale /Jr different from /jf. 

Suppressing the functional dependences, the evolution equations for the singlet parton distri- 
butions read 



Pqq Pqg 



gg 



(2.2) 



where C?) stands for the Mellin convolution in the momentum variable, 



^ l'^aiy)b(-) . (2.3) 



y ^y^ 

The quark-quark splitting function iqq in Eq. ( 12.21) can be expressed as 

^qq = Pns+MPm + Pi^^ = ^nt+^ps • (2-4) 

Here P^^ is the non-singlet splitting function which we have recently computed up to the third order 
in Ref. The 0{a^) quantities and P^^ are the flavour independent ('sea') contributions 
to the quark-quark and quark- antiquark splitting functions Pq^q^ and / q,qj, respectively. The non- 
singlet contribution dominates Eq. (12.41 ) at large x, where the 'pure singlet' term P^^ is very small. 
At small X, on the other hand, the latter contribution takes over as xP^^ does not vanish for x — > 0, 
unlike -^^^ns • The gluon-quark and quark-gluon entries in Eq. (12.21 ) are given by 

-^qg ~ '^fPqig ' -^gq ~ -^gq,- (2.5) 

in terms of the flavour-independent splitting functions iq.g = iq.g and Pgq. = Pg^. . With the excep- 
tion of the al part of Pqg, neither of the quantities xP^g, xPgq and xPgg vanishes for x^O. 

Our calculation is performed in Mellin-A/^ space, i.e., we compute the singlet anomalous di- 
mensions Yab(^7 OCs) which are related to the splitting functions by a Mellin transformation, 

Yab(A^,as) = - / dxx^-^Pab{x,as) . (2.6) 
Jo 

The additional relative sign is the standard convention. Note that in the older literature an additional 
factor of two is often included in Eq. (12.61) . 
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The calculation follows the approach of Refs. lE51E51I^Hm . The optical theorem and the 
operator product expansion are employed to compute the Mellin moments of (partly fictitious, 
see below) deep-inelastic structure functions. Since the moment variable is now an analytical 
parameter, we cannot apply the techniques of Refs. Il25ll26ll . where the integrals were solved using 
the Mincer program 11411 l42l . The introduction of new techniques was therefore necessary, and 
various aspects of those have already been discussed in Refs. Il40ll43ll44ll45ll38l . A salient feature 
of our method is, however, that we can check our extensive manipulations at almost any stage by 
falling back on a Mincer evaluation of fixed low-integer moments. Note also that we will obtain 
the three-loop coefficient functions in DIS as well, once the present calculation is supplemented 
by a second Lorentz projection required to disentangle the structure functions and F^ ll46l . 

The complete set of NNLO singlet anomalous dimensions can be extracted from the third-order 
amplitudes of the forward Compton processes 

parton (P) + probe (Q) — > parton (P) + probe (Q) , (2.7) 

where the probes are the photon (y) and a fictitious classical scalar (|) coupling directly only to the 
gluon field via (SfG^y,G^ . The inclusion of the latter, required for obtaining also the anomalous 
dimensions ygq and ygg to the desired accuracy, leads to a substantial increase of the number of 
diagrams as shown in Table 1. Among the partons in Eq. (12.71) we also include an external ghost h. 
This is done in order to allow us to take the sum over external gluon spins by contracting with 
—gfjw instead of the full physical expression which would, due to the presence of extra powers of 
P, lead to a complication of our task. For similar reasons we do not keep the gauge dependence in 
our all-// computations, but check its cancellation only for a few fixed values of N. 



process 


tree 


1-loop 


2-loop 


3 -loop 


qy qy 


1 


3 


25 


359 


gy gy 




2 


17 


345 


hy gy 






2 


56 


q^ 




1 


23 


696 


g^ g^ 


1 


8 


218 


6378 


— i> 




1 


33 


1184 


sum 


2 


15 


318 


9018 



Table 1 : The number of diagrams for the amplitudes employed for the calculation of the three-loop 
singlet anomalous dimensions. The roles of the ghost h and the scalar (|) are discussed in the text. 

The diagrams are generated automatically with the diagram generator Qgraf ll47ll . For all 
symbolic manipulations we use the latest version of Form P^l49ll . The calculation is performed 
in dimensional regularization Il50ll5ni52ll53ll . The renormalization is carried out in the MS-scheme 
Il54ll55l as described in detail in Ref. ll25ll . using the result of Refs. |l56l|57l| for the renormalization 
of the operator CLG^ entering the scalar case. 



3 



3 Results in Mellin space 



Here we present the anomalous dimensions JahiN, tts) in the MS-scheme up to the third order in 
the running coupling constant a,. The N"LO expansion coefficients y^',^^ (A'^) are normalized as 

' as ^ "+!..(„ 

11=0 



(3.1) 



The anomalous dimensions can be expressed in terms of harmonic sums 11^171 15^115^1501 . Recall 
that, following the notation of Ref. ll5Sll . these sums are recursively defined by 



i=l 



and 



S± 



mi,ra2,...,mi. 



;=1 



(3.2) 



(3.3) 



The sum of the absolute values of the indices rrik defines the weight of the harmonic sum. Sums up 
to weight 21— \ occur in the /-loop results written down below. 

In order to arrive at a reasonably compact representation of our results, we employ the abbre- 
viation Sff, = Sff,{N) in what follows, together with the notation 



N±5„-, = 5,^(A^±1) 



N±i5r„ = Sa{N±i) 



(3.4) 



for arguments shifted by ±1 or a larger integer i. In this notation the well-known one-loop (LO) 
singlet anomalous dimensions dO read 



,(0), 







Y^'^^(A^) = 2n^(N_+4N+-2N+2-3)Si 
Yi?(A^) = 2C^(2N_2-4N--N+ + 3)5i 
yfgW = Q(4(N_2-2N_-2N+ + N+2 + 3)5i-y)+^n^ 
The corresponding second-order (NLO) quantities |I51|6l[T0l[TT|| are given by 



(3.5) 



y'pJiN) = 4C^n^f^(N_2-N-)5i-(N 



(N--N+) 



25i-t-52 + 253 



+ 



56, 8, 
yS, + -52 



+ (1-N_ 



85i-4S2 



(3.6) 



25i+52 + 253 



-(N+-N+2) 



Si 



44 



+ 45i,i + yS2 +(1-N+) 275i+ 451,1-752-253 - 2(N_ +4N+ - 2N+2 - 3) 



218 



1,1 



+ 4C^n^(2(N+-N+2) 



+ (N--N. 



751 --5, 



55i+25i 1-252 + ^3 



43 7 

ySl+45l,l--S2 



+ 2(N- + 4N+-2N+2-3) 



1 



Sl,l,l — 51,2 — 52,1 + 2*^3 



(3.7) 



Y^q^A^) = 4QC^(2(2N_2-4N_-N++3) 



Si. 



1,1 



-S1-2 — S17 — S2. 



-1351,1-752-253 
+ (N+-N+2 



22 

^5i-y5i,i 



44 8 
y5i + -52 



+ (N_2-2N_+N 

+ 4C^«^((N_2-2N_+N 



+ 4(N_-N+) 



4 20 

351,1 -y5i 



+ (1-N+) 

7 
9 

-(1-N 



25i 

5i+352 + 53 



45i 



'25i,i 
•353 



+ 4C/ (2N_2 - 4N_ - N+ + 3) 



-5i+252 + 253 



35i,i — 25i,i 1 



(l-N+)[5i-25i,i + -52 



(N--N+) 



(3.8) 



2 16 23 14 2 

y£)(A^) = 4Qn^(---5i--(N_2 + N+2)5i + y(N_ + N+)5i + -(N_-N+)52) 



.3 3 

+ 4C/{ 25-3 - I - ^5i + 253 - (N_2 - 2N_ - 2N+ + N+2 + 3) 



45i,_2 + 45i 2 + 452,1 



+ ^ (N+ - N+2)52 - 4(N_ - 3N+ + N+2 + 1) 



352-53 



+ — (N_+N+)5i + y(N^ 



- N+)52) + 4Qn^ (^^ + ^(N_2 - 13N- - N+ - 5N+2 + 18)5i + (3N- - 5N+ + 2)52 
-2(N_-N+)53) . (3.9) 



The pure-singlet contribution (12.41) to the three-loop (NNLO) anomalous dimension yqq{N) is 
Yg^(A^) = 16QC^«^(^(4N_2-N--N+ + 4N+2-6)[35iC3 + 5i,-2,i-5i,i,-2 + 5i,i,i,i 
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571 6761 3^ 52^ 56^ 20^ 



8^12 

1^351,-3 + 25i,3 + ^"^1,1,1 + 3*^2,1,1 



+ (Nh 
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19 2 

2Si — 5ij + —52 — -52,1 



(1-N_ 



25] — 5ij +52,1 + 2^5*2 



+ (N_+N+-2) 



77 



63 



37 



1 



L54^i-54^u-2752 + -5i,i,i 



+ -(N--Nh 



5 29 

^52,1 -52,1,1 + -^53 



■253,1 —54 



+ 16C^2^^((N+ 



+ 452,2- ^52,1,1 +35i,2-y 54 
+ 5l,l,l,l-^5i,3] +(N_2 + N+2 



-N+2 

1 



16 13 163 85 28 22 

T^^'' + 54^^ - -n^' - 12^^'' + -9^^' - -J^' 



+ (1-N- 



-5'2,i,i +553,1 



^ (4N_2 - N- - N+ + 4N+2 - 6) [35i ^3 - 5i,i,i - 5i,i,2 
55^ 523^ 23^ 1 55^ 46^ 523^ 

TTr'^l-TTTTT'^l,!-— "^1,2 - y5i + y5i,2 + — 5i,i 

5„ .„ 181 



12 108 ' 9 

298 121 2707 497 63 

- — 5i,2 + ^5i,i - —52 - -52,1 + ^5i,i,i + 552,2 + ■ 

47 971 275 755 5 

"^1,2 — TTTTT"^!,! + TTTT'^l ~ -^^^2 — —"^1,1,1 +652^3 —^'2,3 



-53—54 



+ (N_-N+) 



72 



12 



108 ' 216 

32 3 
+ 1752,1 +252,1,1,1 —252,1,2 — 352,1,1 +252,2 — ^"^s — 253,i — 53,1,1 +453,2 — -54 + 654,1 



■455 



(3.10) 



(2) 

The non-singlet part of Tcv({N) can be found in Eq. (3.7) of Ref. [38|. The third-order results for 
the off-diagonal anomalous dimensions JqgiN) and JgqiN) in Eq. (12.21) are given by 
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— 45*2, _3 +4^2,-2 g"'^2,l —352,1,1 —252,1,1,1 +52,1,2+ 252,2 + 552,2,1 -4^2,3+653,1,1 

+ 1 153 - 453,-2 + y 53,1 - 653,2 - y 54 - 754,1 



+ (N--N+) [^51 + ^51^3-^51,2 



35^ 103^ _ 7^ 13^ 9^ 7^ 1 

+ 35i,-3 - y5i,_2 g~ I'l ~45i,i,_2 - g5l,l,l - y5i, 1,1,1 + -5i,i,2+ -5i,2,l - -52,1,1,1 

9 1 27 3 o 87 

- 25l,3 + ^52 -352^3 +752,-2+ -g-52,1 + ^52,1,1 +52,1,2-52,2,1 +352,3 - Y^53 -53,1,1 



13 27 7 

y 53,1 + 253,2 + y 54 + -54,1 - 355 



71 19 

- 1 15i,_3 + y 5i,-2 + 125i,i,_2 - -g 



+ (1-N 
13 



1759^ 63^^ 17^ 
1 /5i,i - yy5i - y5iQ3 + y5i, 1,1,1 



13 3 13 . ^ 

■5i,i,i — y5i,i,2 + y5i,2 — 25i,2,i + y5i,3 — 352^3 



409 59 1 3 565 

^y 52 — 452,-3 —52,-2 + y52,i — ;^52,i,i + -52, 1,1,1 —352,1,2 + 352,2,1 —552,3 + ^y^a 



17 

— 853,-2 + y 53,1 + 353,1,1 — 653,2 



103 21 

y— 54 — y54,l + 1 155 



(3.11) 



and 



yg^iV) = 16QC^n^((2N_2-4N_-N+ + 3) -^Si - 251^3 + i5i,-3 + -r^5i,_2 - ^5i,3 



967 



41 



1 



144 



3 18 3 

2c 251 4^ 13^ 5^ 5^ 10^ 5^ 151 

351,-2,1 + y^5i,i - -5i,i,_2 - y 5i,i,i + -5i,i,i,i - -5i,i,2 + y5i,2 - -5i,2,i - —52 



1 10 5 1 

^52,-2 + 7r52,l - 752,1,1 + -52,2 

J y D i 



+ (N--N+) 



331 28 11 

— 51-452,-2 + y5i,-2--5i,i,i 



4 2 53 733 4 22 
+ 353,1 - ^52,1 + -5i,i - —52 + -52,1,1 - y 53^ 



+ (1-N+) 



10 1 1 

Ly^2'-2+ 1^^2,1-^51,1 



17^ 137^ 5^ 1 565^ ^ 35^ 2^ 

- T^i'-2 - T^^i + 6^1,2 + 451,1,1 + —52 - 52,1,1 + y^53 - -53,1 



:(N--N+2) 



53 



131 



■4N 



4 

•N+ + 3) 



352,1 H — ^5i +5i,_2 — -^5i,i — 5i,i,i + -^52 
163 „ 3„ 3„ 6503 



125 



32-51 -25i,-4-25i,-3+ ^32 



-(N_-1)54) + 16QC/((2N_2 



5i,i —55i, -2,-2 — 35i. -2,1 — 45i,i,i,i,i 



+ 5i,-2 + 25i,_2,i,i -95i,i^3 -45i, 1,-3 +35i, 1,-2 + 25i, 1,-2,1 +55i, 1,3 + 65i, 1,1,-2 + 5i, 1,2,1 
35 2 1 

+ 35i,i,i,2 + y5i,i,i,i + -5i,i,i - — 5i,i,2 



191 ^ 41 ^ ^ 

y^5i,2 — 35i,2,-2 — Y^5i,2,i +45i,3 — 452,1 



+ 25i,2,i,i - ^5i,4 - ^52,1,1 + 252,1,1,1 +52,1,2 + 352,2+52,2,1 - 252,3] + (N- - N+2) 



652. 



173, 26, 2, 335, 7, 28, 8,1 ,,r, 

+ ^^1^1 - T^i'i^i - S^i'i'i'i - + 25l - 252,1,1 - y 53 + -54J - 6(N_ - 1) [52,-3 



8 



- 252,1,-2 + 3S2C3] + (N- - N+) [36S1C3 - + 125i,_3 - 365i,_2 - ^^1,1 + 4^3,2 

- 165i,3 - 245i,i,_2 - -^Si^i,i + -5"!, 1,1,1 - — '^i,2 + 25i,2,i + 432 "^z + 3^*^2,-2 + 



+ Y^'^2,1 - ^"^2,1,1 - 352,1,1,1 + ^-^2,2 - 2^2,2,1 



119 47 

—53 - 1453,-2 + y53,i -753,1,1 +455 



+ 1052,3 

2525 



+ (1-N_ 



2005^ 117^^ 39^ 315^ 

Si 2~ ~ ^ '^■^ ~^ "J 1,-2 — J 1,-2,1 +35i,i,i,i — 254,1 



55 197 11 53 13 

+ lU^''' +4^^M,-2 - Y2^1,l,l - 35i,i,2 + — 5i,2 - y5i,2,l + y5i,3 + y53,l,l - 452,2 
^^52 - 3752,-2 + 1353,-2 + 2'^2,1,1 + 2-^2,1,1,1 - ^"^3,1 + 352,2,1 - 1252,3 + ^^"^3 



3 ^ 57 

+ 2*^2,1 -653,2 - y54 



+ \6C/Cp f (2N_2 - 4N- - N+ + 3) 



138305 
L 2592 



"^i — 25i,-2,i,i 



11 49 ^ 109 3 ^ 3379 

^'^1,-4 + -^5'i.-3 +5i.-2.-2 — 105i, 1,-2,1 + -rz-5i,-2 — -5i,-2,l +25i 2,2 — ^rr^'^i,! 

Z Iz z zlo 

+ 85i -3,1 +35i,i^3 + 125i.i -3 + — 5i, 1,-2 + 25i, 1,1, 1,1 + — 5i,i,i — 65i,i,i,-2 — ;r5'i, 1,1,1 

z z4 D 

55 71 55 

— 45i,i,i,2 + Y2'^i.i'2 — 45i, 1,2,1 +25i,i,3 + 24'^i'2 + 55i,2,-2 + "^1.2,1 — 45i,2,i,i +65i,2,2 

11 3 395 11 

+ ^5'i,3 +45i 3,1 — -5i,4 — 3—52 — 752,-3 — -^^'2,-2 + 452,-2,1 +252, 1,-2 — 252,i,i,i 

Z Z 34 D 

+ y52,l,l + 352,1,2 - ^^2,2 + 352,2,1 - 352,3+453,1,1 - 453,2] + (N- - 1) [652^3 - 852,-2,1 



+ (N_-N+) 



57595^ 31 143^ 25^ 

^ 5i — I251C3 — 5i —3 ^52 —2 H — ;r5i —2 1 

L 1296 ibJ g 1, J ^ z, z ^ I, 2,i 



689^ 50^ 

^Sl,l + y5u,-2 

31 



11 11 229 113 2200 

+ Yg'^l.l.l ^ y^l, 1,1,1 + ^'^1,2+^51,3 - y^52 -352,-3 - 1253,2 + 95i,-2 + y52,l 

13 37 25 463 

- 1852.1,-2 + ^.52,1,1 +452,1,1,1 - -^^2,2 - ^^2,3 -3153-953,-2 - —53,1 +453,1,1 +54 
D i z Iz 



- ^"^4,1 -855 



37^ 2^ 
+ y5i,i,i + -5i, 1,1,1 - 

77^ 

— TTr53 — 653 1 + ^54 
18 ' 3 



+ (N--N+2 
145 



4 2105 8 109 4 

i^S] —9 1 : — 5i — — 5i —3 — 105i —9 — - — 5i i — — 5i i _2 

3 " 81 3 ' 27 ' 3 ' ' 



4 584 104 8 14 

J^5i,2 - -5i,3 - —52 - 452,-2 - -^'^2,1 + 3-^2,1,1 - y ^2,2 



+ (1-Nh 



39 ^ . 29843 ^ 17 ^ 145 ^ 29 ^ 



25^ 57^ 13^ 5^ _ 97^ _ 41 7417^ 

- y^i,-2 ^ y^i,i,-2 - ^^"^1,1,1 + 4*^1, 1,1.1 +45^1,1,2 - ^5^1,2 + 451,2,1 - y^i,3 + 

1 92 53 9 1 

+ 2*^2,-3 + y 5*2,-2 — ^^"^2,1 + 1552,1,-2 — ^^2,1,1 —352,1,1,1 +552,2 + -54,1+3853 + 853,2 



41 9 92 25 31 

+ y5'2,3 + 2*^1,-3 + ^"^1,1 ~253,i,i + ^"^4 + y^s 



+ 16C^n/(-(l-N 



2^1 — 5i,i 

2 c 371 
3^^'2-432^^ 

35 1 1 1057 16 8 1 2 181 2 

y5i,-2 - -5i,i - ^-^i.i.i + + y^^'"^ ~ 9*^^'^ ^ 3"^^'^^^ ~ 3"^^'^ ^ ~\2^^ ~ 3"^^'^ 



^(2N_2-4N--N+ + 3) 
6 



2^1 + -5i,i — 5i,i,i 



+ 16C^2„^((N_-N_ 



9 



3^4 + 455 
1625 



r 1 31 95 1 1 

+ (2N_2 - 4N_ - N+ + 3) [25iC3 - 351,2,1 - j^Su^2 + + -Si,3 + 3^1,2 



144"^i- 6^1-1,1,1-3^1,1,2 



7 c 83^ 2^ 

YQg5i,i+-5i,i,i + -52,-2 



i(N--N+2) 



■S1-2-S3 



+ (1-N- 



15137 49 



107 19 5 

-^5i,i + — 5i,i,i - -Si,2 - 1052-452,-2 
io li b 



1 155 

2'^2,l,l +52,2 - ^^53+53,1 +54-655 



+ \6C/ (2N-2 - 4N- - N+ + 3) 



65i. 



^2.-2 



47 7 47 

■ttSi — 5i.-4 — -5i 2 + 65i 3 — — 5i,i +65i,i^3 +45i 1 3 — 65i 1 2 — 35i 1 2 — 35i 1 3 
Id z Id 

23 9 7 

-^•^l,!,! — 2'^1'1'1'1 +25i,i,i,i,i +5i,i,i,2 + 35i, 1,2,1 + -5i,2 + 25i,2,-2+25i,2,l,l — 25i,2,2 



+ 2(N-- 1) 652^3 -452,1,-2 + 853,-2 + (N_ 



287 
'32 



5i -2451^3 +5i,i,i,i 



- 125i,-3 + 365i,-2 + ^5i,i + 165i,i,-2 + ^5i,i,i + ^5i,2 - 25i,2,i +95i,3 - 453,i - 552,3 

91 41 35 

+ 353 1.1 — —52 + 852,-3 — 3052,-2 — ^"^2,1 +52,1,1 —52,1,1,1 + 252,2,1 — —53 — 54 + 354,1 
Id 4 o 



r 749 141 433 17 

+ (1 -N+) 3951C3 - — 5i +205i,-3 - ^^1,-2 - -r^5i,i +65i,i,i - — 5i,i,i,i 
L d4 z Id 4 



-255 
— 305i 1,-2 — 5i 1.2 



19 3 57 9 37 

-jSi2 + 2'5l,2,l - -jSi3 + 2152 - 1052,-3+3552,-2 - 2'^3,i,i + ^-^4 



19 9 3 11 485 27 9 

+ ^^2,1 + 4-^2,1,1 + 2-^2,1,1,1 + 352,2 - 352,2,1 + y 52,3 - ^^3 + —53,1 - -54,1 



. (3.12) 



Finally the three-loop gluon-gluon anomalous dimension reads 



Yg(iV) = 16QC^n^(^|^ + (N-2-2N--2N++N+2 + 3)[45iC3 ^ 

2 4 521 16 1 4 4 17 8 

- 3^1,-3 + 351,-2,1 - y^5u - -5i,i,-2 + -5i,i,i - -5i,i,i,i + -5i,i,2 - -5i,2 - -5i,3 

86 4 2 2 4 

70 31 7 7 55 1 

- Y^i,i,-2 + ^5i,i,i - -5i,i,i,i + -5i,i,2 - -5i,3j + (N- - N- 



/241 



15331 44 

5l — -;r'^l.-2 



+ (N- +N+ - 2) [l75iC3 + y5i,-3 - -5i,-2,i 

133 221 

L7^5i,2-^5i,-2 



673 4948 49 119 
~ '54^^^^ ^ ~Sr^^ ~ ~ ^^^^^^ ~ ^^^'^^ ^ ^^^^'^^ ^ T^^'^ ^ 16^2,1,-2 + 654,1 

7 251 10 29 

— 75'2.i,i +252.1.1,1 —252,1,2 — 52,2 + 752,3 + -r7r5'3 — ^5'3,i — 5'3,i.i +453.2 — —54 + 855 

D Iz i D 



-8(N_- 1)53,-2 + (N--N+2) 



127^ 511 97^ 103^ 

53 5] — 65] —2 5] 1 — 35] 2 + 253 I 57 

18 12 ■ 12 ■ ' ' 27 



8„ 16^ 2^ 

3^2,-2 -y52,i- -52,2 



+ (1-N+) 



1807^ 604^ 5311 52^ 1667 
-324^' + ^^1,-2 + ^5i,i - -5i,2 - ^ 



S2 



68 53 7 19 67 9 33 

^52,-2 - y^2,i - 3'^2,l,l + "^"^2,2 + Y^53 + -53,1 + ^54 - 2055 



324 

2c 44^ 4^ 521 16^ 1 4^ 4^ 8^ 

+ 3-^1,-3 + y^i,-2 - 3*^1,-2,1 + Y^-^i.i + y^i'i -2 - Q-^i.i.i + 3'^i, 1,1,1 - 3'^i.i.2 + 3*^1,3 



10 



17^ 86^ 4^ 



17^ 

13 29 

+ 54^1,1 - + (^-2 - 2N_ - 2N+ + N+2 + 3) 
-252,1+53 

+ (N_2 - 2N- - 2N+ + N+2 + 3) 

2^ 205 
3 

o r 31 4 
— o5i^3 — — 5i -3 + -5i, 
D i 



^^2,1 - ^52,1,1 + ^52,2) + 16Qn/ (ii(l - N+)52 - ^5i 
I - 2N+ +N+2 + 3) [^5i - i^5u] - ^(N_ - N+) [52 



--(N_-N+) 52 



,rl204 , 2 19 11 

-j^Si - 45i^3 - 3^1,-3 + ^^1,-2 + 25li,-2 + y 5i,2 



81 

- 3^1,-2,1 + J^S,,i - -52 - -52,-2 + y52,lj 

31 4 2441 4 25 1 

y Si,_3 + 35,-2,1 + ^^1,1+ 9Sl,l,-2 + ^5i,2 + -5i,3_ 

59 71 2 3 64 5 

+ 652^3 + 352,-3 - "5"'^2,-2 - J2^^'^ ~ ""^2,1,-2 " 35'2,2 " 2*^2,3 - -gSi + 553,-2 + ^2*^3,1 

- 254 - ^S4,l] + (N- - N+2) [^52,-2 - ^52 + y 53 - ^53,l] + (1 - N+) [^52 + 55 



+ (N-+N+-2) 
25 



305 ^ 1405 ^ 

"'^1,-2 TTTT'Jl 



18 



648 



+ (N-N,) -S, 



109 



25^ 8^ 473^ _ 25^ 31 i 10^ 1 5443^ . 

+ ^'J2 —2 ~ 77'J2 1 T^'JS ~ 453 —2 ^53 1 H — —54 — —Si 3 TT— "J 1 + 25i C3 

2 



3 
2 

3^1-3 - y^i -2 -r 3^1-2,1 - Yos^i.i ~ ^^1,1-2 - y^i:2 -r 3 

y ^2,1 - y 53 - ^53,1) + 16Q3 _ _ 2N+ + N+2 + 3) 

88 85 



^ 37 2 205 13 2 151 2 

+ :j^l,-3 - y^l -2 + 3"^! -2,1 - YoS"^!.! ~ 25i,i,-2 - y5i,2 + 3'^-2,-2 + ^^2 + 3'52,-2 

\ »x r -5i-165i,-4 



9 " ' 3 

648 

- 115i,-2,i +45i,-2,2- ^5i,i +245i,i,-3 

67 55 

+ 115'i,i,-2 — 165i, 1,-2,1 + 85i,i,3 — y5i,2 + 85i,2,-2 + 85i,2,2 + y^i.B + 85i,3,i — 85i,4 

11 67 22 

-3 - y S2,-2 + 852,-2,1 - y S2,l + 452,l,-2 + 852,1,2 + y 52,2 + 852,2,1 

61 



88 85 

+ -r'^1,-3 + 165l,-3,l + ^5i,-2 + 45i,-2,-2 
3 D 



395^ IJ 

52 - 1452-3 

27 ' ^ 3 



+ (N- + N+ -2) [l45i,-2,i - ^5i - y5i,-3 - y5i,-2 



- 1052,3 + 853,1,1 -853,2 

80 109 1 , 

- 2^^51,1 + 145i,i,-2 - — 5i,2 + 45i,3j + (N- 

- 852,1,-2 + ^^2,2 - 1052,3 + ^^3 - 2053,-2 + ^53,1 - 1653,2 - 2154 - 2654,1 

3 3d 3 

/XT XT ^\o 9533 77 ^ ^ 44 1517 
+ (N_ - N+2) 852,-3 — Jos' '^2 - y 5'2,-2 - 852,-2,1 - 852,i,-2 - y^2,2 f8~'^3 - 855 



r473 

-N+) ;77T'52- 1252,-3 +552,-2 -252,1 
Lzlo 



+ 853,-2 3 



-^53,1+453,2 + 4454+ 1654^ 
109 



+ (1-N+) 



r8533„ 103^ 1579^ 



98 
-y54 



79 0/ 

- :^ + 45-5 - 85-4,1 + y 5-3 - 45-3,-2 - 25-3,2 - 45-2,-3 ■ ^ 

88 523 22 



67 



67 



413 



32 

11 16619 

^5-2,-2+45-2,-2,1 +45-2,1,-2 



-^^1,2+108^1,1 



33 781 

- 1 15l,l,-2 - y 5i,3 + —52 - 452,-3 



11 

+ y52,-2 + 452,-2, 



67 

1 — 7r52,i +452,1,-2 I ^ 
9 o 



11 



+ ^53,1 



11 



+ y 53 - 453,-2 - 253,2 - 854,1 + 455) + l6Cpn/ { (N_2 - 2N_ - 2N+ + N+2 + 3) 



77^ 16^ 2^ 
~8T^^'^27^^'^~9^^'^'^ 



+ -(N_+N+-2) 



11 2^ 16^ 

y^ + ^5i,i,i--5i,i 



211 139 11 

— IS""^^'^ + y^2 + 5'2,i +'^2,1,1 —252,2 — 253,1 +54 



11 



64 



58 



1 



(N- - N+2) [25i -5i,i + -52 + -52,1 - ^53] + (1 -N+) [-5i + -5i,i + -53 



10 1 

y52 + 3^2,1 



1 



+ l6C/nf (^3 (N_2 - 2N_ - 2N+ + N+2 + 3) [-5i,2 + 2'5i,3 - ^1,1,1 



31 



-'5'i,-3 + 25i,i,_2 + Y^5'i,i +5i, 1,1,1 



11 
16 



"^l — 5'l,l,2 



2 

+ (N-+N+-2)[^5i,3-95iC3 



16^ 67^ 23^ 7^ 7^ 32^ 

yJi,-3 + yJi,-2- Y^"^!,!,! + 3 "^1,1, 1,1 ~ 3'J1,1,2 + y5i,i,_2 



+ (N_ 



■N^ 



254 1 —255 



~ 3*^1, 1 + -^'^2 + 652^3 +452,-3 - -^^2-2 " 35'2,1 - 852,1,-2+ 35'2,1,1 +252,1,2 



- 252,1,1,1 - y52,2 - 352,3 - y 53 - 453, 1 +.^3,1 J + + y5i,2] + (N- - N+2) 



85 
12^ 



'1,1 



163^ oc 9^ 8^ 4^ 4^ 4^ 14^ 2^ 
+ ^y^l ~ 35i,2 - -52 + 3-^2,-2 - 352,1 + 3-^2,1,1 - 352,2 + y 53 - -54 



+ (1-N- 



454 



^y5i,i — 85i,2 + y 52 + 852,-2+ y 52,1 +5'2,1,1 —352,2 — jy5'3 — 5'3J + y^i 



1 11 4 31 8 4 4 4 

+ 8(N_ - 1)53,-2 -J^ + J^Si + -5i,-3 - -5i,i - -5i,i,-2 + -5i,i,i - -5i,i,i,i + -5i,i,2 



-3^1,2-3 



^1,3) 



(3.13) 



Eqs. ( 13.101) - (13.131) represent new results of this article, with the only exception of the C^nJ part 
of Eq. (13.131 ) which has been obtained by Bennett and Gracey in Ref. iFTl . Our results agree with 
the even moments N = 2, ... ,12 computed before |E51E^ using the Mincer program IHTIH^ . 

The results (13.51) - (13.131) are assembled, after inserting the QCD values Cp = 4/3 and Q = 3 
for the colour factors, in Figs. 1 and 2 for four active flavours and a typical value tts = 0.2 for the 
strong coupling constant. The NNLO corrections are markedly smaller than the NLO contributions 
under these circumstances. At A'^ > 2 they amount to less than 2% and 1% for the large diagonal 
quantities Yqq and Ygg, respectively, while for the much smaller off-diagonal anomalous dimensions 
Yqg and Ygq values of up to 6% and 4% are reached. The relative NNLO corrections are very large 
at A'^ > 2 for Yps, which is however completely negligible in this region of N. 

For — i> 00 the off-diagonal n-loop anomalous dimensions vanish like — ln^"^^A^, while the 
diagonal quantities behave as 



(A^) = Al {InN + Y.) - - C + O ( - 



A^ 



A^ 



where Ye is the Euler-Mascheroni constant. The leading large-A'^ coefficients An of Yqq have been 
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Figure 1: The perturbative expansion of the diagonal anomalous dimensions JqqiN) and JggiN) 
for four flavours at = 0.2. The pure-singlet (ps) contribution to Yqq is shown separately. 
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Figure 2: As Fig. 1, but for the off-diagonal anomalous dimensions Yqg(A^) and ygq{N). 
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specified up to n = 3 in Eq. (3.1 1) of Ref. As expected, the constants are related to those 
results by 

Al = . (3.15) 

The coefficients C| in Eq. (13.141) can be expressed in terms of the Af by 

cf = , C| = 4CAAf = (Af)2 , cf = SCaAI = 2A\aI . (3.16) 

This result is completely analogous to the corresponding relation for Cn in Eq. (3. 12) of Ref. ll38ll . 
Finally the A'^-independent contributions Bn can be read off directly from the 5(1 — x) terms in 
Eqs. Clll, dOOh and (RTTSl) below. 



4 Results in x-space 



in) 

The N"LO singlet splitting functions P^^ [x) in 



i'ab(as,^)= L (4.1) 



;?=0 



are obtained from the A'^-space results of the previous section by an inverse Mellin transformation 
which expresses these functions in terms of harmonic polylogarithms ll63l l64l l65l . This trans- 
formation can be performed by a completely algebraic procedure l65ll based on the fact that 
harmonic sums occur as coefficients of the Taylor expansion of harmonic polylogarithms. 



Our notation for the harmonic polylogarithms //mi,...,m„,('^), mj = 0, ±1 follows Ref. 165| to 
which the reader is referred for a detailed discussion. For completeness, we recall the basic defini- 
tions: The lowest- weight (w = 1) functions H,n{x) are given by 

Ho{x) = \nx , H±i{x) = ^ ln(l ^jc) . (4.2) 

The higher- weight (w > 2) functions are recursively defined as 

{^In'^jc, if mi,...,m^ = 0,...,0 

(4.3) 
dzfmi{z)H^2,-,mM) ^ else 

with 

/oW = -, /±iW = T^- (4.4) 
For chains of indices zero we again employ the abbreviated notation 

^0,...,0,±1,0,...,0,±1,...(^) = H±(^m+l),±{n+l),...{^) ■ (4.5) 



m n 
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Corresponding to the maximal weight 21 — 1 of the harmonic sums in section 3, the /-loop splitting 
functions involve harmonic polylogarithms up to weight 21—2. Hence our three-loop results cannot 
be expressed in terms of standard polylogarithms which are sufficiently general only for w < 3. 

For completeness we recall the one- and two-loop non-singlet splitting functions 





and 



/201 



2nfPqg{x) 

2CpPgq{x) 

Q(4pggW + y5(l -jc)) - ^/i^5(l -;c) 



3 " 9 



+ {l+x] 



5Ho - 2H 



0,0 



(4.6) 



(4.7) 



4C4«/ (y ~ -2 + 25x- 2pqg(-;c)H_i,o - 2pqg(x)Hi j +x~ 



+4(1 -jc) 



-^2 



+ 4x' 



Ho,o — 2Ho+xHi 
5 



Ho + Hon + ■ 



-4^2^-6Ho,o + 9Hoj +4Cpnf[2p^g{x) 

29 



44 218 



Hio + Hij -I-H2 



+ 2(1 Ho + Ho,o-2xHi + 



y-Ho,o-^Ho) (4.8) 



P^!,\x) = 4QC^(-+277gq(x 
-7Ho + 2Ho,o-2Hix+(l+;c) 



Hi,o + Hi,i+H2-^Hi 



44 1 



-Ho 

3 9 



2Ho.o — 5Ho + 



37 



^H, 



+ 4Cp-(pgq{x) 



9 J 

3Hi-2Hi,i 



+ 4^2- 

'2pgq(-x)H_i^oj -4Cpnf(^^x 



+ {l+x) 



Ho,o-^ + ^Ho 



+ I--H0 + 2H1X 



-3Ho,o 

(4.9) 



11 



2 Ho + 8Ho,o - 2 



+(l+x) 
44 9 

-yX Ho + 2/?gg(jc) 



4^)(x) = 4Q«^(l-x-^pgg(x)-^(i-x2j_2^1^^)H^_2g^l__^^^j^4^^2(^^ 

27 



13 /I 



+ 2pgg(-x) 



H 



0,0 



'2H_i,o — ^2 



2, 

3' 
67/1 



9 Vjc 



x^) -12Ho 



^ - C2 + Ho,o + 2Hi ,0 + 2H2] + 5( 1 - x) + 3^3 



+ 4C^n/2Ho 



21 10 o 



4 - 5Ho - 2H 



10,0 



^5(1 -x)) . 



(4.10) 



Here and in Eqs. (14.121) - (14.151) we suppress the argument x of the polylogarithms and use 



Pqg(x) = \ —2x + 2x 

Pgq{x) = 2x^^ —2+x 

Pgg{x) = (1 —x)^^ +x^^ —2 + x — x^ 



(4.11) 
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Divergences for x ^ \ are understood in the sense of +-distributions. 



The third-order pure-singlet contribution to the quark-quark splitting function (12.41) . corre- 
sponding to the anomalous dimension (13.101) . is given by 



,(2), 



Pps'{x) = l6Cj^Cpnf(-{-+x 



4,1 



13 14 1 , 

yH_i,o - yHo + 2^-1^2 - H_L-i,o - 2H_i,o,o 



-H 



-1,2 



+ li-.-^^) f^^2 + H2,1 +9^3 + ^Hi,o 



-(- 

3 X L 3 
3Hi,o,o + 2Hij,o + 2Hi ij 



1 Hi H H? -|- Hi C? Hi 1 

216 72 ^ 3 ^ 6 ' 



+ {l-x) 



182 158 397 13 

-^Hi + + ^g-Ho,o - ^H_2,o + 3Ho,o,o,o 



13 119 

+^Hi,o + 3xHi,o + H_3,o + H_2C2 + 2H_2.-i.o + 3H_2.o,o + xHo,oC2 + ;rHi ^2 - 7Hi,o,o 
6 2 ' 2 4 

3tt tt tt 1 .1 Nr7„. 31^ 91^ 71^ 113. 826^^ 
--Hi,i +Hi,i,o + Hi,LiJ +{l+x) [y^HoC2 + y ^3 + Y^H2 + — H3 + —^2 - ^Ho 

5 16 31 17 117, 2 . 5 , 

+;rH2,o + -rH-1,0 + 6xH_i,o + — Ho,o.o - ^H2.i + — ^2 + 9HoC3 + :^ii-i^2 + 2H2,i,o 
z i D zU z 



1 7 1 

+ 2 1 ,0,0 - 2H_i,2 + H2C2 - 2 H2,o,o + H_ 1 ^_ 1 ,0 + 2H2, 1, 1 + H3, 1 - -H4 



+ 5H_2,o + H2,i 



1 32 29 235 511 97 33 

+Ho,o,0,0 - 2 ^2 V 4H_3,o + 4H0C3 - y Ho,o - j^Ho " ^^^2 - ^ - y^Hi + — H2 - H3 

11 113 7 r83 243 51 1 Q7 4 

--H0C2 - y ^3 - ^H2,o - lOHo.o.o + [y Ho,o - y^Ho + 10^2 + ^ + y Hi - -H2 

-4^3 - H0C2 +H3 +H2,o - 6H_2,o] ) + 16C^n/ (— Hq - 2 - H2 C2 + [H2 - ^2 + 3 



-yHo 



2/1 2^ 



+ {l-x) 



1 7 35 185 

-H,,-gH,+xH, + -H„ + — 



+5(1 +^) [jHz - ^C, + C3 + H2,i - 2H3 + 2H„C2 + ^Ho.o + Ho,o,o] ) + 16Cf -n/(^H, 



25 583 101 73 73 

-y Ho,o - Ho,o,o + ^Ho - ^ + y ^2 - y H2 + H3 - 5H2,o - H2,i - H0C2 



55 
12 



85 22 109 13 28 28 16^ . , 16„ , ,4 26 

"Y^ni - y no,0 ^ + —L,2 - y H2 - y HoC,2 + y^s +4hl2,o + 3H2,l - y 1,3 



22 

+yHo,o,o 



4/1 2, 



23 523 55 1 



+(l-x) 



1 7 2743 53 251 5, 5 8 

L2Hi,o,o + Y^Hi,i - yy Ho - -Ho,o - —Hi - -^2 + 4H2 - -Hi,o + 3xHi,o 



1669 5 , 37, 2 

^TTT + xHoAO +4H2,i +7H2,o + 10x^3 - 777^2 
zlo z 1(J 

(4.12) 



+3HoC2-3H3-Hi,i,o-Hij,iJ +(l+x) 

— 7Ho^3 + 6Ho,oC2 — 4Ho,o,o,o + H2,o,o — 2H2,i,o — 2H2,i,i — 4H3^o — H3 j — 6H4 
Due to Eqs. (13.111) and (13.121) the three-loop gluon-quark and quark- gluon splitting functions read 



d(2)( 

qg ^ 



-39 



15, 



leC^CpnApqgix) — Hi^3 -4Hi,i,i +3H2,o,o - -rHi,2 + tHi,i,o + 3H2,i,o 



4 ' 4 
2 49„ 3, 



1 



+H0C3 - 2H2,i,i + 4H2C2 - ^HoC2 - ^Ho,o + y - C2' - y H2 - 2Hi,o,o,o - 3Hi,o,o 



16 



-^Hi,o - yHi,i - —Hi + -H2,o + 2Hi^2 + yHo,o,o + - ^ + 

+6R2A +3Hi _2,o +9Hi^oC2 + 6Hi^i^2 + Hi,ho,o + 3Hij^i^o — 4Hi^i^i^i — 3Hi^i^2 — 6Hi^2,i 

49 1 rl7 5 5 9 5 3 

-6Hi,3 + —^2 +7'qg(--«) yH_i^3 - -H_i,_i,o - 2^-1.2 - 2^-1,0 + 2^-^fi + 2^-i,0fi 

-2H3,i - 2H4 - 6H_2,2 + 6H_2,-i,o - 6H_2A0 + 2Ho,oC2 + 9H_2C2 + 3H_k_2,o - 2H_i,2,i 
-6H_i,_i,_i,o + 6H_K_i,o,o + 6H_i,_i,2 + 9H_i,oC2 - 9H_l_i^2 - 2H_i,2,o - yH_i,o,o,o 



1,3 



,1 2.155 23^^ 4^^ 

+^^-^HT2-^^^+y^^^«-3^^'^'«. 



— H| on H| H Hi 1 

3 ' ' 108 9 ' 



2 

-3H,,,, 



+ (1 -j) [6H2.i,o+3H2.i,i - 5hi.i,i -7H2.o,o-2Hi,2 + 39HoC3 -4H2C2 - ^Cs 



154 ^ 899 121 ^ 2 607 5 , 65 29 13 

+Hi,i,o + — H0C2 + ^Ho.o + — C2' + ^H2 - -Hi^2 + -^Hi,o,o - y^Hi,o - — Hij 

1189 67 949, 67 142 215 3989 

- - - - -36^' - y - + ^ - + 

+(1 +x) [h_i,o,o - 10H_2C2 + 6H_2,o,o + 2Ho,oC2 - 9H_i,_i,o - 7H_i,2 - 9H_2,o - 2H3,i 

37 5 1 

-4H_2,-i,o - 4H4 - 4H3,o - 4Ho,o,o,o + yH-1,0 + ^ " 4H_2,o,o + 2Ho,oC2 

9 11 19 9 

+H2C2 — 3Hij^o + 2Ho,o,o,o + H_3^o — 9H2j,o — 21^2,1,1 + -^^1,1,1 + yH2,o,o + 2^^'^ 

91 5 5 9 39 473 1853 
— 2^o^3 + 8H_2C2 + 2^-1-1,0 + 2^-1,2 + 2^-^^o + ^^-2,0 - ^Ho^2 48~ 1^0,0 

217, 59, 2 169 13 , 2 167 191 1283 185 

75^^ 170, 85^^ 425^^ 7693 3659^^ ^ 
H H2 H C2 H Hn H H3 H 1 Hn - 2x 



xH2,2 + 4H3n-4H_2,2 



1 



229 



11 



+ 16Qny (^-pqg(x)|^Hi,2-Hi^2-Hi,o,o-Hi,i,o-Hi,i,i-— Ho + -Ho,o+ 2j 

53,, 17,, ^ lU 139] 1 , „, 53 J 2 2 

~T^Ho + — Ho,o-^3 + -r^^2-77^ +:TPqg(-^)H_i,o,0- T^l---^ 



+ x 



1 



H2 



18 

7 7 
--xHi — Hq.o + T-^Hi4 
o z 



18 



108 



162';c 



:(1-^) 



L6 
6Ho,o,o 



7 7 19 5 5 

+ -x{l +x)H_i,o + -Ho - —Hi +Ho,o,o + ^Hij + -H_i,o 



+ 16C/n^(pqg W [3Hi,3 + ^Hi,o,o - yH2,i + ^^2' - ^Hu,o + ^H3 - yHiC3 

5 63 23 155, 25 2537 867 23 

-^H2, - y H,o - ^H,2 - + ^H2 - ^Ho + — - y H_i,o,o + 3H4 - Hi,i, 

383 25 3 7 ^1 103 5 2561 

+ ^Hi^i - yH-2,0 - g 4)2 - ^HiC,2 - itlo,04)2 - ^^^04)2 + y^^i + 2^''^''^''^ ~T2r ' 

+Hi 1 1 — 2H2,o,o — 3Hi 2,0 — 5Hi^oC2 + 3Ho,o,o — Hi4^2 — Hi j,o,o — 4Hi r j n + 2Hi j r i 



-2Hij,2 — 2Hi^2,o 
727 



+ Pqg 



36 



H_i,_i^2 - 2H_L2 - 6H_i,_i,o + Hi,i,i + 2H_2C2 - H_2A0 
+— H-1,0 - H_i^2 - 2H_2,2 - ^H_iC3 - H-1,-2,0 + 2H_i,_i,o,o + 2H_i,_i,2 - ^H_i,o,o,o 



17 



+6H_i,_i,_i,o - 2H_i,3 + 2H_i,2,i 

161 2351 

H 1 — 



8 3 ^161 2351 

-3H-,„,„ + -H,,„+6C3 + 3^H,- — 

1 n 



1 .2^^2H ,32 4 IC 



-^^1,1 



-2H 
3 

"2 



3 5 31 

2243 265„ 33^^ 31„ 23^^ 497^ 29^^ ^ 143^^ 

^ + ^H_,o,o + y H2,o,o + 19H2a + ^2^1^ + yH-2,0 - ^^2 + ^Hi^^ - 

IItt 19^^^ 1223„ 43„ 3011^ 1 ,^ ,r 



L3 

2,1 o,r26 28 

1,2 + Hi^2 + H_i^2 + +Hi,i,i] +{\-x) [l5Ho,o,o,o -5H2C2 - yCs + 
' "^H +^^H ^^^2 29 113 18691 

■n2,0 + -j^ni,0-yj-4)2 - yni,0,0 ^"^2 + 



72 



-Ho 



11 19 1223 43 

--^Hi 1.1 - — Ho^2 + ^;:;-Hi - — Ho.o.o — ::rr~^o.o 
6 12 72 



+ {l+x) 8H2,i,o-4H_i,2 



O IZ /Z JO J L ' ' 

+7H_i _i^o — -^Hi^a — 5H_2^2— llH_2,o,o+ 2^-1,0 + yH_i^2 + 8H3a — 10H_2_] 
+5H2C2 + 4H2,i,i - H_3,o + 36H0C3 - 5H2C2] + 2H„i,2 + 6H_i,_i,o - 6H2,i,o - 3H2,i,i 



25 13 27 11 13 17 

-1 lHo,o,o,o - 5H3,i + — Hi^a + yH_2Q2 + yH_2,o,o + yH_3^o + y H2Q2 - yHi^o.o 

17 3 1, 11 79 67 263, 2 

+ 13H_2,-l,o - J2^W - 4H4 - 4Ho,oC2 + Hi,2 + y Hi, 1,0 + —Uifi + yHi,o + — ^2^ 

,119 967 305 o4Hr+H r ^^^^^H ^^^^ 38H ^^H 

217 7 79 4 17 17 31 

- ~ 2^-^'^ ~^ j2 3^^ 12^^'^'^ 12^^^ 18^^ 3Ho,o,o 



3 
79 

'4 
145 



__H2oo- — 

4'' 24 z /z J iz iz 10 

7„ 11^ 739 163^^ 7„ 

— Hi \ Ho H Ho + 2Ho o o o 

36 96 24 24 ' ' ' ' 



1553,^ \ 2/7„ 11„ 739 163 

■H3 H Ho + 16Cpn/ -Ho H Hi \ , 

^ 24 ' / V6 ' ' 36 96 24 24 

5,, 1 , , r 91 35 



5 5 5, 1 , , 

^H2--Hl,o+-C2 + ^PqgW 



12 

5 5 5 

-5H,,.-5H.--H, 

7 1 77 1 
-C3-2H100 + -H1 +— (- 

7 
9 



77,1 2n 



yv zt zt 

91 35 22 

H2a + y - y Ho - y Ho,0 + Hi,i,i + 6Ho,0,0 

6463 ^„ 16^ 7 ^ 

- yHo,o,0 + g-^Hi4 



7 8 7 

+ -xH2 + -JcHi,o- -JC^2 



— Hi 4Hn n n 

,12 432 " ' 3 



,[3475^^ 103^^ 

— 2H_3,o — 7Hi ^3 + 5H2,2 + 6H3,o + 6H3 j + H2, 1 ,0 + 4H2,o,o + 3H2, 1 + 2H2,i, 1 + 2 ^2,0 
61 61. 87„ 11„ 61„ 17„ . 5„ 5„ 



[o + ^Ho,o ^ + l6Cp-nf(^pqg{x) 7Hi,3 + 7H4 



z 

ui 61, 87 11 61 17 , 5 5 19, 

+^H2 - y C2 + y Hi + y Hi,2 + y Hi,i + y Hi,o " 7Ho,oC2 + 2Hi,o,o + 2^1,1,0 - y ^3 

y Ho + y ^2^ + 3Hi,i,i - 5H2C2 - 7H0C3 



81 11 11 , 7 , 15 87 11 

+32 + y - y^«^^ - 2^^^^ + y + + y - - 

o,o-5Hij^2 + 4Hi,i,o,o + Hi,ij,o + 2Hij,ij+5Hij,2 



32 ' 2 ^ 2 " 2 ' ■ 2 ' 
+ 1 lHo,o - 2Hi,_2,o - 7Hi,oC2 + 3Hi,oao 

+6Hi,2,o + 6Hi,2,i +4pqg(-x) Ho,o,o,o - H-2,0 + H-i -1,0 - H_2,o,o + ^H_i _2,o - ^ 



^H-in 



5 11^ 1 1 0/ 

-^H_i,o,o- 2^-3,0+ 2^"i ^2 + H-i, -1,0,0 - ^H_i,o,o,oJ +2(1 -jc) |^H2,i,o-H2,o,o-H2,2 

9 43 49 13 
-H3,i - 2H3,o - 2H_i^2 +Hi,2 - Hi,o,o - Hi,i,o + H2^2 - ^2" + y ^2 + "g"^2 + y^i.i 



18 



33 5 7 21 479 1 1 1 1 3 



17 19 239 405 

+ 2H0C3 - 2^4 + 1^1^2- ^Ho,o,o- ^^Ho^o- ^2" ^0 



+ 8(l+;c) 



H 



H 



-1,0,0 



3 9 

-Ho,o,o,o + ^H_2,o - 4^-1,0 



113 71 35 11 33 7 7 5 

+4xH_2,o,o g~^2 - -^^2 - "^Hi - ^Hi^2 - "g"^i>i ~ 2^^'^^ ~ 2^^'°^^ ~ -j^l-Q-Q 

-2^1,1,0 - 2^3 - - J J - -H3 - -H2,i - 2^2,1,1 + ^Ho^2 + H2^2 + 2H0C3 

9 2^7 49 391 401 

+7C2 + + -Hi ^2 + -rHcco + ~rF^o,o + -tf^o — H2,o,o — H2,i,o + H2,2 + H3,i 
3 z z 4 Id Id 

1 



+2H3,o + 6H_i^2 + 2H2,o + 2H_2C2 + 4H_2,-i,o 



(4.13) 



and 



d(2), 

gq ^ 



\6CjiCpni^ 



25 131 125 

—Hi - + 3^2 - H_i.o - 3H2 + Hi.i + —Ho - Hq.o 

D 4 D 



5 r 967 251 39 2 1 4 

+ g/'gq W [Hi,2 + H2,i + + ^Hi - -Hi,i - 3^3 - 3H0C2 - 3H1C2 - 3 Hi,o + Hi,i,o 



2 2 
"^Hi^o,o + Hij,i + ^^2,0 



2 r 7 41 151 1 

+ ^Pmi-^) [2H_iC2 + 4^2 + Y^H_i,o - —Ho + -H_2,o 



+^H2 + 2H-i,_i,o - H_i,o,o - H_i,2] + ^(1 --^) [H-2,0 + 2^3 - H3] + (1 +x) 



179, 
108 



Hi 



5^ 25 5 167 1 4 

+ n^2 + 7rH-i,o — TzHi.i — TF^O.O — o H2,i — -H0Q2 
y y JD io i i 

1„ 227„ 35„ „ 2„ ^ 10 



193 1 1 

^ + 4«i + 9«-i'0 + 4H2 



^ 2 rll 

4^^'^ T8~^° ~ 12^^'° ^ H2,i - -Ho^2 + ^H_2,o + 3^3 + 2H3 + -Ho,o,o +^ -^^2 



523 19 271 5 

rl? H rln ill n 

144 36 108 6 ' 



173, 



26, 



+ 16QC^^(x^ [- + ^Hi -2^3 - 3Hi,i,i - -Hi,i 



335 28 8 1 r3 163 

-6H2 + 2H2,i + 6^2 + ^Ho - y Ho,o - 3Ho,o,oJ + /'gq W [-Hi^3 + ^ - 5^2 + ^C3 

6503 2 35 ^ ^ ^ o r 41 

H — r^Hi + ^Hi^i + — Hij^i +4H2 + 2^2,1 +4Hi^o,o + 2H2,o,o — H2Q2 + ^2^1.2 + ^^2,2 



432 

191 3 59 5 

■^Hi^o + 3H2,o - 2H2,i J - -H_i^2 - Y^^i^2 + 5Hi _2,o + Hi,oC2 + ■:^i,q,q,q - 2Hi j^2 



+^Hi 4,0 + 5Hi 4,0,0 - 3Hi,i ,1,0 - 4Hi ,1,1,1 - Hi,i,2 - 2Hi,2,i +H2,i,o 



+ Pzq{-x) 



H-i,o 



3 27 2 11 3 

+H_i,oC2 + 2^-1,0,0 + ^^^2 - 3H_i,_i,o - — H_i^3 - 3H_i,_2,o - - 3H_i,2 

+5H_i,_i^2 - 4H_i,_i,o,o - 2H_i,_i,2 + 6H_i,_i,_i,o + 2H_i,2,il + (1 -x) [H2C2 - H2,2 



,23 7061 4631 38 4433 21 

12 '° ^ I32" ~ T44^ '° ~ Y '"■'^ " ^-^fi ~ ^^3,0 - ^^^1 ~ ^H2,o,o - yHi,o,o 



1^2' - ^Hi,2 + yHi^2 - 4H0C3] + (1 +^) 



— hl3 - hl_2,0 - -^no'=2 - 2 36~ 



19 



655 151 185 1 95 29 



171 



H-Lo — 12H_i,o,o + 7H_i^2 



5 3 

+ 1 6H_ 1 ^_ 1 ^0 + 3 ^2,0 + 2 ^2, 1 , 1 + 4Ho,o,o,o 



4 

35H_2,0 - + ^^^0,0 - y Ho,o,o 

13 , 13 15 2005 157^ , 1291 55 

-2H3,o - yHoC2 - 13H_3,o - yH3,i + yH3 - + — C2 + 8^3 + + -Hi,i 

3 1 27 11 3 5 

+-H2 + -H2,i + yH_i,o - y Hi,o,o - 8H2A0 - 4^2^ + ^Hi,! - H2,2 + -Hi^i + 8H_i,_i,o 



+4H2,o + 2^2,1,1 - H_i^2 + 7H2C2 + 6H_2C2 + 12H_2,-i,o - 6H_2A0 +^ 



3Hi 



1,1 



' Ho,oC2 



9 35 
+ 2^-1,0,0 — ^Hi^o + 2H4 + 3Hi j,o + H_i^2 



+ l6C/Cpix^ 



2 2105 77 



16, 14 2 , 14 104 4 37 

-6H3 + y ^3 - 10H_i,o - yH2,o " " y Ho,o,0 + y-H2 - -Hi,o,0 + yHi,i 



4 104 8 145 4 2 

H — H_i _i n C? H? I H Hi n H — H_r 9 H — Hi 1 1 

3 ' ' 9^3 ' 18 -3 '3 ' ' 27 



Hi + -H_i^o,o + 6H0C2 



584 ■ 
+4H_2,o + y-Ho 



, , r7^^ . 138305 1„ 13^^ . ^„ 11„ 
+ Pgq W [2^1^^ + -2592' ~ 3^'° ^ xH_i^2 + 2H2,i,i + yHi,o,o 

43 109 17 71 11 21 3 

+4H3^i - yHi,ij - ^^2 - ^^2,1 - 24^1=0 - "^^-2,0 - y ^3 + 2^1,0,0,0 - Hi -2,0 

395 55 55 

^^-tUq — 2Hi^oC2 — Hij^2 — T^Hi^i^o + 2Hij^o,o + 4Hi^i^i^o + 2Hijjj +4Hij^2 — Y^^i>2 



54 ^ " 12 

+6Hi,2,o +4Hi,2,i +4Hi,3 + 3H2,i,o + 3H2,2 
109„ „ . 17. 2 1 



■23, 



L 2 

\ — r;rH-i,o + Ho^3 + -^^2^ + 7Hi^2 + 2H2^2 — ^ttHij — — H_i,_i,o — 4H3,o — 3H2,o,o 



+ Pgqi-^) + 5H_2C2 + 2H_2,-1,0 

65„ 19, 



12 



5 6 ^ " 24 ' 2 

3 3379 49 11 

-7H_2A0 - :^H_i,2 + — r^Hi - 4H_2,2 - ^H_i,o,o - ^H_i,o,o,o - 13H_i,_i^2 - 8H_i,3 
z lib D z 

— 6H_i _i _i^o + 12H_i _i,o,o + 10H_i _i,2 + 10H_i,o^2 + 5H_i _2,o — 2H_i,2,o — 2H_i^2,i 



+^HoC2l+(l-^) 

97 10 245 1 
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Finally the Mellin inversion of Eq. ( I3.13I ) yields the NNLO gluon-gluon splitting function 
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The large-x behaviour of the gluon-gluon splitting function igg {x) is given by 
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The constants A| and C| have been specified in Eqs. (13.151) and (13.161) . respectively, while the 
coefficients of 5(1 —jc) are explicit in Eq. (14.151) . The corresponding limit of the gluon-quark and 
quark-gluon splitting functions is 
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It is worthwhile to notice that all the coefficients in Eqs. (14.181) and ( 14.191 ) except Df^ vanish 
for the choice 

Ca = ric = Cf =nf (4.20) 

of the colour factors leading to a A'^= 1 supersymmetric theory. This is part of a general structure. 
The combination 



As (x) = P^^> (x) + P^^ (x) - {x) - pS {x) 



(4.21) 



24 



of the (n+l)-loop MS splitting functions is found to be much simpler than the functions P^^\x) 
themselves. In fact, after transforming to the dimensional reduction (DR) scheme respecting 
the supersymmetry, Ag(jc) vanishes for both the unpolarized ||9| and polarized (spin-dependent) 
ll66l l67l l68ll two-loop splitting functions. We are not (yet) in a position to present this scheme 
transformation at the third order. However, we do obtain the above-mentioned simplification within 
the MS scheme; especially all harmonic polylogarithms of weight four cancel in the combination 
(I4.21I) for choice (14.201) of the colour factors. We plan to return to this issue in a later publication. 



We now return to the end-point behaviour. At small x the three-loop splitting functions read 

(4.22) 
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The coefficients of the \/x terms of Pqq (which are, of course, entirely due the pure-singlet con- 
tribution given in Eq. (14.121) ) are given by 
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or, after inserting Ca = 3 and Cf = 4/3 and the numerical values of ^2 and ^3, 
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The corresponding results for the gluon-quark splitting function (14.131) are 
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The coefficients E\ in Eqs. (14.231 ) and (14.251) agree with those obtained by Catani and Hautmann 
in Ref. lETI from the small-x resummation. 



The small-x coefficients of the quark-gluon splitting function (14.141) are given by 
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Finally the corresponding coefficients of the three-loop gluon-gluon splitting function (14.151) read 
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= 2675.85 + 157.269/7^ 
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(4.30) 



The coefficient Ef^ is identical to the result obtained from the next-to-leading logarithmic BFKL 
equation by Fadin and Lipatov in Ref. lESll after transformation to the MS scheme (as given, e.g., 
in Eq. (4.7) of Ref. ll^ ). Numerically the simple relation CpEf^ = C^^f is broken by less than 
2% in the part and less than 0.5% for the complete coefficients at /jy^ = 3, . . . , 6. 



The three-loop splitting functions (14.121) - (14.151) are shown in Figs. 3 - 6 for /Zy = 4 together 
with the approximate expressions inferred in Ref. [ 37 1 from the fixed-A'^ results of Ref"s. Il25ll26ll and 
the small-x limits of Refs. Il27ll28l . Also displayed are the respective leading small-x contributions 
Ef'x'^lnx. Notice that all splitting functions have been multiplied by x for display purposes. 

(2) 

With the exception of Pgq , where no small-x 'anchor' was available, our exact results comply 
with the error bands of Ref. OTll for the full range of x shown in the figures. Hence it is reasonable 
to expect that an extension of the results of Refs. lE51 to the next order, using a future four- 
loop generalization of the Mincer program iHTl H^ . would, together with small-jc constraints, 
facilitate relevant estimates of P^^, {x). We expect that such an extension, while still a formidable 
task, will be performed much earlier than the fourth-order version of the present calculation. 
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Figure 3: The three-loop pure- singlet splitting function (14.121) for four flavours, multiplied by x 
for display purposes. Also shown is the uncertainty band derived in Ref. OTIl using the lowest six 
even-integer moments ll25ll26ll and the leading small-x term ll27l . The latter contribution is shown 
separately on the right-hand-side (dotted line) for x < 0.01. 
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Figure 4: As Fig. 3, but for the third-order gluon-quark splitting function specified in Eq. (14.131) . 
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Figure 5: As Fig. 3, but for the three-loop quark-gluon splitting function (14.141) . Note that in this 
case the leading small-x contribution was unknown before the present calculation. 
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Figure 6: As Fig. 3, but for the third-order gluon-gluon splitting function specified in Eq. (14.151) . 
This diagonal quantity has been additionally multiplied by ( 1 — jc) . The leading small-x term (again 
shown by the dotted line on the right-hand-side) has been first obtained in Ref. 
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Figure 7: The convolution of the three-loop gluon-gluon splitting function (14.151) with schematic 
'steep' (left) and 'flat' (right) gluon distributions. Also shown the results obtained by instead using 
only the leading (Ei) and the leading and next-to-leading {E[ + E2) small-jc terms in Eq. (14.221) . 
and by supplementing the latter by a constant term restoring the correct second moment. 



As also illustrated in Figs. 3-6, the leading small-.)c terms ~ x \nx alone do not provide 
good approximations of the full results (I4.12I) - (I4.15I) at experimentally relevant small values of x. 
At X = 10^^, for example, they exceed the exact values of Pj^\x) by factors between 1.6 and 
2.0 for njr = 4. Good small-x approximations of these quantities are obtained by including all 

x^^ contributions as specified in Eq. (14.221 ) - (14.301) . However this does not apply, as obvious 

(2) (2) 
from Fig. 7, to the convolution [igg' ®^](x) by which igg'^ enters the evolution equations (12.21) . 

Even if the two terms explicit in Eq. (14.221) are (non-uniquely) supplemented by an jc-independent 

contribution restoring the correct second moment, even the sign of the convolution remains wrong 

down to X ~ 10^^ for the simplified, but not unrealistic gluon distribution ~ x^^^^{l —x)^. 

As our exact expressions (14.121) - (14.151) for the the functions P^^^ (x) are neither particularly 
short nor especially simple, we also provide compact approximate representations built up, besides 

(2) 

powers of x, only from the +-distribution (for Pgg (x)) and the end-point logarithms 

©0 = 1/(1-^)+, Li=ln(l-;c), Lo = lnjc. (4.31) 

(2) 

Inserting the numerical values of the QCD colour factors, ips in Eq. ( 14.121 ) can be represented by 
P^^^ (jc) ^ \nf(- 5.926 - 9.75 1 - 72. 1 1 Li + 177.4 + 392.9 jc - 101 .4jc2 - 57.04 LqLi 
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-661.6Lo+131.4L^-400/9L^+ 160/27 L^-506.0;c"'-3584/27x"% 

+ nj 1. 778 + 5.944 Li + 100.1 - 125.2x + 49.26x2 - 12.59x3 - 1.889LoLi 

+ 61.75 Lo+ 17. 89L^ + 32/274 + 256/81x"i)}(l-jc) . (4.32) 

Correspondingly the off-diagonal quantities (14.131) and (14.141) can be parametrized by 

pf^{x) ^ n^(^100/27L|-70/9L^-120.5L? + 104.42Li+2522-3316x + 2126x2 

+ LqLi ( 1 823 - 25.22 Lq) - 252.5 xL^ + 424.9 Lq + 88 1 . 5 - 44/3 
+ 536/27 - 1268.3x^1 - 896/3x"%) 

+ nj (^20/27Lf + 200/27 L?-5.496Li- 252.0+ 158.0;c+145.4;c2 

- 139.28^3 -LoLi (53.09 + 8O.6I6L0) - 98.07 xL^ + 11.70xL^ 

- 254.0 Lo - 90.80L^ - 376/27 - 16/9 + 1 1 12/243x^1 ) (4.33) 

and 

pf^{x) ^ +400/81Lt + 2200/27Lf + 606.3L? + 2193Li -4307 + 489.3^+1452x2 
+ 146.0x3 -447.3 L^Li - 972.9xL^ + 4033 Lq - 1794 L^ + 1568/9 
-4288/81L^ + 6163.1x"^ + 1189.3x"% 
+ n^iy- 400/81 - 68.069 - 296.7 Li - 183.8 + 33.35x - 277.9x2 

+ 108.6xL^-49.68LoLi + 174.8Lo + 20.39 + 704/81 
+ 128/27L^- 46.41x^1 +71.082x"1Lo) 

+ nj (^96/27Lf (x^^- 1 + 1/2 x)+320/27Li(x"^- 1+4/5 x) 

-64/27 (x"i - 1 -2x)) , (4.34) 

where the nj part is exact. Finally the gluon-gluon splitting function (14.151) can be approximated by 

P^^^{x) ^ +2643.521 ©0 + 4425.8945(1 -x)+3589Li -20852 +3968X- 3363x2 
+ 4848 x^ + LqLi (7305 + 8757 Lq) + 274.4 Lq - 747 1 + 72 - 144 
+ 14214x"^+2675.8x"% 
+ nfi^ -412.172 ©0-528.7235(1 -x) -320Li -350.2 + 755.7x-713.8x2 

+ 559.3x3 +L0L1 (26.15 -808.7Lo) + 1541 Lo + 491.3L^ + 832/9 
+ 512/27L^+ 182.96x^1 + 157.27x"%) 

+ n^(^- 16/9^0 + 6.46305(1 -x) - 13.878 + 153.4x- 187.7x^ + 52.75x3 
-LoLi(115.6-85.25x + 63. 23 Lo) - 3. 422 Lo + 9.680 - 32/27 

- 680/243 x"^) . (4.35) 
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The coefficients of 1 jx, (Inx) jx, \v? x and In^x are exact in Eqs. (14.321) - (14.351) . up to a truncation 
of the irrational numbers. The same holds for the coefficients of ln-^(l — x) and ln'^(l — x) in 
Eqs. (14.331) and (14.341) . and those of ©o and ln(l — x) in Eq. (14.351) . The remaining terms (except, 
or course, for the 5( 1 — ;c) parts in Eq. (14.351) ) have been obtained by fits to the exact results (14.121) 
- (14.151) at 10"^ <x< 1-10"^ which we evaluated using the Fortran code of Ref. |l69|. Smaller 
values of x are not needed, as all Xjx terms are exact. Except for values of x very close to zeros 
of P^^{x), the parametrizations (14.321) - (14.351) deviate from the exact expressions by less than one 
part in a thousand, which should be amply sufficient for foreseeable numerical applications. 

Finally the coefficients of 5(1 —x) in Eq. (14.351) have been slightly adjusted from their exact 
values using the lowest integer moments. This is a somewhat tricky point, so let us briefly elab- 
orate on it. For Pqq and /gg the low moments, and partly also the convolutions with the parton 
distributions, involve large cancellations between the integrals over the (fitted) regular parts and 
the 5(1 —x) contributions. The second moment of the Hy^-independent part of /gg, for example, 
vanishes due to the momentum sum rule (recall that iqg has no contribution, cf. Eq. (12.51) ). 
while the respective third-order coefficient of 5(1 —x) is as large as 4 ■ 10^. A maximal accuracy 
of the parametrization (14.351) . and of the convolutions with the gluon distribution, is thus achieved 
by 'fitting' this coefficient to the second moment. For the case under consideration this actually 
leads to a very small adjustment of about 0.01%. 

One important approach to implementing higher-order results into the numerical evolution of 
the parton distributions and the analysis of general hard processes is the moment-space technique 
iTTUl 171117^17^ IT?! , which requires the analytic continuation of the anomalous dimensions (12.61 ) 
to certain complex values of A'^. Also these complex-A/^ moments can be readily obtained to a 
perfectly sufficient accuracy using Eqs. (14.321) - (14.351) together with the corresponding non-singlet 
results in Eqs. (4.22) - (4.24) of Ref. The Mellin transform of these parametrizations involve 
only simple harmonic sums S,n>o{N) of which the analytic continuations in terms of logarithmic 
derivatives of Euler's F-function are well known. The reader is referred to Refs. ||5^ 17511 for a 
more mathematical approach to the analytic continuations. 



5 Numerical implications 

We are now ready to illustrate the numerical effect of our new three-loop splitting functions P^-^ {x) 
on the evolution (12.21) of the singlet-quark and gluon distributions qs{x^ijj) and g{x,pj). For all 
figures we choose a reference scale = — 30 GeV^ - a scale relevant, for example, for deep- 
inelastic scattering both at fixed-target experiments and the ep collider HERA - and employ the 
sufficiently realistic model distributions 

xq,{x,n^) = 0.6x-°-\l-x)^-^l+5.0x^-^) 

xg{x,ni) = 1.6x"0-3(l-x)4-^ (1-0.6x^-3) (5.1) 
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irrespective of the order of the expansion. This order-independence does not hold for actual data- 
fitted parton distributions like those in Refs. Il33ll34ll . but here it facilitates direct comparisons of 
the various contributions to the scale derivatives / = dlnf/dln/jj for / = t/s, g. For the same 
reason we employ an order-independent value for the strong coupling constant, 

asiMo) = 0.2 , (5.2) 

corresponding to a fairly standard value at the Z mass, as(Mj) ~ 0. 1 16, beyond the leading order. 
Finally our default for the number of effectively massless flavours is «^ = 4. 

The respective scale derivatives of the singlet-quark and gluon distributions are graphically 
displayed in Figs. 8 and 9 over a wide range of x. Numerical values can be found for four charac- 
teristic x-values in Tables 2 and 3, where we also show the dependence on rijr and the break-up into 
the quark- and gluon-initiated contributions. As these two terms can occur with different signs, 
and since the LO and NLO results partly display a somewhat anomalous behaviour (see below), 
the picture is much less clear-cut here than in the non-singlet sector discussed in Ref. ll38ll . 
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Figure 8: The perturbative expansion of the scale derivative = dlnq^/dln/jj of the singlet 
quark density at /jj = /Jq for the initial conditions specified in Eqs. (15.11) and (15.21) . 

The scale derivative of the quark distribution (Fig. 8 and Table 2) is dominated at large x 
(small x) by the iqq ® (Pqg®g) contributions. The former (latter) is actually negligible for very 
small (large) values of x. The NNLO corrections are small at large x with respect to both the total 
derivative and the NLO contributions. At small-jc all NLO contributions are very large (or the LO 
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terms are abnormally small, recall that xP^ and x/qg"* vanish for x — > 0). Consequently the total 
NNLO corrections, while reaching 10% at x = 10^^, remain smaller than the NLO results by a 
factor of eight or more over the full x-range. 
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Figure 9: As Fig. 8, but for the gluon density. The spikes close to jc = 0. 1 in the right parts of both 
figures are due to zeros of the LO and NLO predictions and do not represent large corrections. 

The situation is rather different for the evolution of the gluon density (Fig. 9 and Table 3). The 
contribution from Pgg ® g dominates for all x (except for extremely large values not considered 
here), but the igq ® part is nowhere negligible. Already the NLO corrections are small especially 
at small x and furthermore the g- and ^j-initiated terms cancel each other to some extent. Thus the 
ratio r2/ri of the relative NNLO and NLO corrections is rather large at small values of x, despite 
the NNLO contribution amounting to only 3% for x as low as 10^^. 

We now turn to the stability of the perturbative expansions in Figs. 8 and 9 under variations of 
the renormalization scale /j^. For /j^ ^ the expansion of the splitting functions in Eq. (14.11) is, 
using the abbreviation = a^/ {4-n) , replaced by 



2/..2^ 



,(1) 



ab 



(5.3) 



+ 



..'a.?) (f<?-{p>''r+2Po/'.';'}ing+p?/'ri„^g ) + 



A/; A*; 

where represent the MS expansion coefficients of the (3-function of QCD 117611771178117911 . 
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= 3 , tts = 0.4 and modified input (see Table 2) 
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Table 2: The LO, NLO and NNLO logarithmic derivatives of the singlet quark distribution at four 
representative values of x, together with the ratios r„ = N"L0/N"-1L0 — 1 for the default input 
parameters specified in the first paragraph of this section and some variations thereof. 
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Table 3: As Table 2, but for the scale derivative (iln^/(iln/j2 of the gluon distribution. 
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In Figs. 10 and 11 the respective consequences of varying /Jr over the rather wide range 
< l-i^ < S/jj are displayed for the logarithmic /j /^-derivatives of the singlet-quark and gluon 
distributions (15.11) at six representative values of x. In both cases the scale dependence is consid- 
erably reduced over the full A:-range by including the third-order corrections. With the exception 
of the smallest x-value considered, x = 10^^ (and of jc = 0.05 in Fig. 11, where the derivative 
is very small anyway), the points of fastest apparent convergence and of minimal /j, -sensitivity, 
df/dfjr = 0, are rather close to the 'natural' choice /jy = /Jf for the renormalization scale. 
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Figure 10: The dependence of the NLO and NNLO predictions for the derivative dlnq^/dlniuj of 
the singlet-quark distribution on the renormalization scale /Ur for six typical values of x. The initial 
conditions are given in Eqs. (15.11) and (15.21) . 



The relative scale uncertainties A^s and Ag of the average derivatives, estimated using the 
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Figure 11: As Fig. 10, but for the derivative g = dlng/dln/jj of the gluon distribution. Notice 
that the scales of the ordinates of the graphs differ within as well as between the two figures. 



conventional interval j/jf < /Jr < l/jf, 

^. ^ max [f{x,Hr = ^Pf--- 2fJf)] - min [fjx.iJr = {pf--- ^Pf)] 
2 1 average [/(jc, = ^A*/ ■ ■ ■ 2a//)] I 

are finally shown in Fig. 12. For the singlet-quark (gluon) distribution, these uncertainty estimates 
amount to 2% (1%) or less at x > 10^^ (4 ■ 10^-^), an improvement by more than a factor of three 
with respect to the corresponding NLO results. Taking into account also the apparent convergence 
of the series in Figs. 6 and 7, it is not unreasonable to expect that the effect of the higher-order 
singlet splitting functions will be about 1% or less for x > 10^^. Larger corrections have to be 
expected at small x. One should also keep in mind that at fourth order also terms with the colour 
structure d'^^'^ dabc/nc — which enter the non-singlet case already at three loops and have a large 
effect dX X < 10^ ll38l — will contribute to the singlet splitting functions. 
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Figure 12: The renormalization scale uncertainty of the NLO and NNLO predictions for the scale 
derivatives of the singlet-quark density (right) and the gluon distribution (left) as estimated by the 
respective quantities A^s and defined in Eq. (15.41) . 



6 Summary 



We have calculated the complete third-order contributions to the splitting functions governing the 
evolution of unpolarized flavour-singlet parton distribution in perturbative QCD. Our calculation is 
performed in Mellin-A/^ space and follows the previous fixed-A'^ computations ll25ll26ll inasmuch as 
we compute the partonic structure functions in deep-inelastic scattering at even N using the optical 
theorem and a dispersion relation as discussed in [25|. Our calculation, however, is not restricted 
to low fixed values of A'^ but provides the complete A/^-dependence from which the x-space splitting 
functions can be obtained by a (by now) standard Mellin inversion. This progress has been made 
possible by an improved understanding of the mathematics of harmonic sums, difference equations 
and harmonic polylogarithms 1158116511401 . and the implementation of corresponding tools, together 
with other new features P^ . in the symbolic manipulation program Form [48 | which we have 
employed to handle the almost prohibitively large intermediate expressions. 

Our results have been presented in both Mellin-A'^ and Bjorken-x space, in the latter case we 
have also provided easy-to-use accurate parametrizations. We agree with all partial results avail- 
able in the literature, in particular we reproduce the lowest six even-integer moments computed 
before ll25ll26l . We also agree with the resummation predictions of Refs. Il27ll28ll for the leading 
small-;c logarithms (Inx) jx of the splitting functions /qq, /qg and /gg, and with the large-n^ result 
ll6Tl for the simple C^n? part of Pgg. Our results respect the supersymmetric relation between all 
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four splitting functions for = = to the extend expected for the MS scheme. At large x 
we verify the expected simple relation between the leading 1/(1 —x)j^ terms of /qq and igg. We 
find that also for the gluon-gluon splitting function the coefficients of the leading integrable term 
ln(l —x) at order n = 2,3 arc proportional to the coefficient of the +-distribution 1/(1 — at 
order n — l,m complete analogy with our surprising findings in the non-singlet case ll38l . 

We have investigated the numerical impact of the three-loop (NNLO) contributions on the 
evolution of the singlet-quark and gluon densities. At x ^ 10~^ the perturbative expansion for 
the scale derivatives / = d\nf{x,/uj)/d\nijj, f = q^, g appears to be very well convergent and 
suggests a residual higher-order uncertainty of about 1% or less at ttg < 0.2. Consequently the 
perturbative evolution can be safely extended to considerably larger values of tts, hence lower 
scales, in this range of x. The situation is much less clear at smaller x. For = 0.2 and realistic 
initial distributions with xqs, xg ~ jc^*^'^ at small x, the NNLO corrections for and g rise towards 
X ^ 0, respectively reaching 13% and —6% at x = 10^^. We stress that the results of the small-jc 
resummation alone cannot help here. For example, not even a qualitatively reliable prediction can 
be expected for the convolution Pgg®g, by which igg enters the evolution equations, even when all 
1/jc terms are included. Besides knowledge of as many of these terms as possible, further progress 
at small x would require at least a four- loop generalization of the fixed-A^ calculations ll25ll26l and 
of the X- space approximations lITTll linking them to the small-jc limits. 

Form files of our results, and Fortran subroutines of our exact and approximate split- 
ting functions can be obtained from the preprint server |http : / /arXiv . orq| by downloading the 
source. Furthermore they are available from the authors upon request. 
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